Local orientational ordering in fluids of spherical molecules 
with dipolar-like anisotropic adhesion. 
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We discuss some interesting physical features stemming from our previous analytical study of a 
simple model of a fluid with dipolar-like interactions of very short range in addition to the usual 
isotropic Baxter potential for adhesive spheres. While the isotropic part is found to rule the global 
structural and thermodynamical equilibrium properties of the fluid, the weaker anisotropic part gives 
■ rise to an interesting short-range local ordering of nearly spherical condensation clusters, containing 

short portions of chains having nose-to-tail parallel alignment which runs antiparallel to adjacent 
f***^ , similar chains. 
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ON ■ Even simple hard sphere fluids display a non-trivial phase diagram, as a function of the packing fraction, which can 
be experimentally probed and theoretically interpreted^. Softening the potential and/or increasing its range, leads to 
a remarkably richer phase diagram which has attracted considerable attention recently (see e.g. Ref. for a recent 
review). Yethiraj and van Blaaderen 3 - have discussed how it is experimentally possible to tune the interactions from 
hard sphere to soft and dipolar ones. More recently, Lu et alA have shown that, contrary to an intuitive expectation, 
gelation of particles with short-range attractions is intimately connected with its equilibrium phase diagram. 

It is widely believed that the addition of a long-range repulsion to a short-range attraction inhibits phase separation, 
by promoting the formation of an equilibrium gel. The same mechanism can be achieved by reducing the probability 
\ of forming a bulk liquid using the concept of limited- valency and/or patchy particles^. This idea has been recently 
explored by a number of authors, who have used the so-called Kern and Frenkel model with circular adhesive patches 
Q ' (of non- vanishing area), or that with short-ranged attractive point-sites on the surface of hard spheres 5 -^^ ' 10 ' 11 . 
O . In spite of their usefulness, the above models share a common shortcoming on the discontinuous dependence of 
the potential on the particle orientations, which makes them very difficult to investigate from a theoretical point of 
view. This drawback is not present in molecular interactions where this dependence is continuous, as for instance in 
dipolar interactions^, a case which is particularly interesting for various reasons. First, because of their widespread 
appearance in colloidal suspensions, such as ferrofluids, which have important practical applications. In addition, 
recent studie d 13 ' 14 ! 15 have shown the existence of a significant influence, in the equilibrium properties of the fluid, of 
chain-like aggregation characteristic of the dipolar interaction, which strongly competes with a stable fluid-fluid phase 
' I . separation. 

Motivated by this features, in this paper we then take an extreme alternative of considering a tail with dipolar- 
like anisotropy combined with a very short-range attraction. The latter is patterned after the well-known Baxter's 
sticky hard sphere (SHS) potential, where attraction occurs only at contact^. Building upon our previous, almost 
fully analytical, study on this model within the Percus-Yevick closure with orientational linearization (PY-OL)^i, we 
discuss here some additional interesting features on the local ordering properties, which were not accounted for in our 
previous work. 

In the same spirit of Baxter's isotropic counterpart^, the model is defined by the following Mayer function^ 
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f(l,2) = fH8(r)+te(l,2)aS(r-a), (1) 

where fns{ r ) = O (r — a) — 1 is its hard sphere (HS) contribution, the Heaviside step function ( Q(x < 0) = 0, 
0(x > 0) = 1 ), and the Dirac delta function <5 (r — a) ensures that the adhesive interaction occurs only at contact (a is 
the HS diameter). The symbol i = (r», fij) (with i = 1,2) denotes both the position r, of the molecular center and the 
orientation f2j which combines the usual polar and azimuthal angles {O^Lpi). Thus we have: (1,2) = (rj.2, f2i, fig) = 
(r,ri2, Ox, O2) = {r, O r , Six, Q2), with rx2 = T2—T1, r = | r 12 1 , and il r being the solid angle associated with rx2 = 1*12/7". 
Moreover, t is the stickiness parameter, equal to (12r)~ in Baxter's original notation^, which measures the strength 
of surface adhesion and increases with decreasing temperature. 

Finally, the angular dependence of the surface adhesion is expressed through the angular factor 

e(l,2) = l + a£>(l > 2), (2) 

including the dipolar function 

D(l, 2) = D(Q r , Six, 2 ) = 3(f • m)(r • u 2 ) - ui • u a (3) 
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Here and in the following, the unit vector represents the orientation fi^ of molecule i, while r coincides with 

r 12 = — r 21- 

The anisotropic function e(l,2), which has the same symmetry as the dipolar interaction, modulates the sticky 
attraction. The requirement e(l, 2) > along with —2 < D(l, 2) < 2 enforce the limits < a < i on the anisotropy 
degree. This range corresponds to the surface interaction always being attractive. In the isotropic case, one has a — 
and e(l,2) = 1. 

As convolutions of Mayer functions generate correlation functions with a more complex angular dependence^ 2 -, it is 
necessary to consider also the angular function 

A(l,2) = Ul -u 2 , (4) 

whose limits of variation are clearly —1 < A(l, 2) < 1. 

We note the difference between the dipolar anisotropic adhesion introduced here and the anisotropy belonging to 
the class of uniform circular 'sticky patches'— £ i 18 i 19 i 20 i 21 i 22 . i n the latter case, the strength of adhesion is uniform, 
independent of the contact point inside an attractive patch, whereas in our model the value of the anisotropic 
correction at D(l, 2) changes with the position of the contact point. Moreover, D(l, 2) can assume both positive and 
negative values, depending on the molecular orientations. Consequently, the strength of adhesion between two particles 
1 and 2 at contact depends - in a continuous way - on the relative orientation of Ui and U2 as well as on the unit 
vector f 12 of the intermolecular distance. The orientations with D(l, 2) > 0, and thus with e(l, 2) = 1 + a D(l, 2) > 1, 
correspond to an attraction stronger than the isotropic one (given by e(l,2) = 1), whereas the configurations with 
D(l, 2) < 0, and thus with e(l, 2) < 1, arc characterized by a weaker attraction, which can even reduce to zero (HS 
limit) in the case of highest anisotropy admissible in the present model, i.e. a = 1/2. 

In particular, we shall focus on a set of parallel and antiparallel configurations with A(l,2) = ui ■ 112 = 1, and 
A(l,2) = —1 respectively. The surface adhesion reaches its maximum value when Ui = 112 = ?i2, which yields 
D(l, 2) = 2 and e(l, 2) = 1 + 2a (head-to-tail parallel configuration). On the contrary, the stickiness is minimum, and 
vanishes for a = 1/2, when Ui = U2 = ?i2, which corresponds to D(l, 2) = —2 and e(l, 2) = 1 — 2a (head-to-head 
or tail-to-tail antiparallel configurations). The intermediate case of orthogonal configuration (112 perpendicular to ui) 
corresponds to -D(l, 2) = 0, which is equivalent to the isotropic SHS interaction. 

Introducing the orientational average (. . .) u = (4ir) J du . . . we note the following results 

(A(l,2)) UiiU2 = (£>(l,2)) UiU2 = 

(5) 

(A(l, 2)^(1, 2)> UiiU2 =0 and (£>*(!, 2))^ = § 

In a previous paper (hereafter referred to as Ref. I)i£, we have analytically solved for this model the Percus-Yevick 
integral equation with an orientational linearization. 

We here recall the main results, referring to Ref. I for details. We start with the molecular Ornstein-Zernike (OZ) 
integral equation for homogeneous fluids 



Ml, 2) = c(l, 2) + pj dr 3 ( c(l, 3) h(3, 2) ) U3 (6) 

where h(l,2) and c(l,2) are the total and direct correlation functions, respectively, and p is the number density. 

Any angle-dependent correlation function F(l,2) could be expanded in a basis of rotational invariants 2 ^, whose 
first few terms are 

F(l,2)=F (r)+F A (r)A(l,2)+F D (r)D(l,2) + --- , (7) 

We stop at the linear terms, assumin g 12 ' 17 that the angular basis {1, A, D} is sufficient for our purposes. 

The PY-OL closure^ is a combination of the PY closure, i.e. c PY = / (1 + 7), with the linear expansion of 7 = h — c 
given by 7° L (1, 2) = 7o(r) + 7a(/")A(1, 2) + 7£>(r)D(l, 2), which also neglects the DA and D 2 terms stemming from 
the product fj. This leads to 

c py-ol (1) 2) = Co(r) + CA(r)A(1; 2 ) + c D (r)D(l, 2) (8) 



co(r) = /Hs(r)[l + 70 (r)] + A a 5 (r - a) 
caM = /hsM7aM + A A 0-5 (r - a) 
cd{t) = /hsMT-dM + k D o-6 (r - a) 



(9) 
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A = t [l + 7o(cr)], A a = t 7a (c), A D = t j D (a) + aA . (10) 
The solution of the OZ equation with the above closure then yields the approximate pair distribution function 

/ Y -° L (1, 2) = 1 + h PY -° h (l, 2) = g (r) + h A (r)A(l, 2) + ^(r)Z)(l, 2) (11) 

' g (r) = ens(r)[l + 7o(r)] + A ad (r - a) 

< h A (r) = e m (r) lA (r) + A A aS (r - a) (12) 
. h D (r) = e H s (r)7u M + A_d(T(5 (r - a) 

where go(r) = 1 + ho(r), and eHs(^) = 1 + /hs(^) is the HS Boltzmann factor. 

The first term in Eqs. ^ corresponds to the well known isotropic Baxter's sticky hard sphere solutionis and the 
OZ equation and this closure constitute a self-contained system. The remaining two have a similar form, but they 
depend in a non-trivial way upon the isotropic term (see Ref. I for details). 

It is instructive to consider the behavior of the <?(12) assuming that r 12 • u x = 1. We focus on a generic reference 
particle 1, with fixed position ri and orientation ui, and consider a particle 2 located along the straight half- line 
which originates from ri and has the same direction as Ui (polar axis). Imagine that 2 has fixed distance r from 
1, but can assume all possible orientations u 2 , which - by axial symmetry - can be described by the single angle 
012 = cos~ 1 (ui • u 2 ). Consequently, 5(1,2) reduces to: g(r, #12) = 9o(r) + [hA(r) + 2/iu(r)] (ui • u 2 ). 

Figure [3a) depicts the behavior of go (r) , which coincides with the reference isotropic part g lsoSHS (r) of the pair 
correlation function, at 77 = 0.4. 

Here, t = gives the HS limiting case, g (r), and we consider increasing values of t, which correspond to increasing 
adhesion or decreasing temperature, i.e. t = 0.1,0.3,0.5 and 0.8. The last t-value yields r = l/(12t) ~ 0.1, which lies 
close to the critical temperature of the isotropic fluid^. 

Two features are noteworthy. First of all, the short-range interactions mainly modify the short-range portions of 
the pair correlation functions. Very pronounced effects are visible in the range a < r < 2a, but significant changes 
are also present all the way out to r — 4cr and beyond, while the phase of the oscillations is clearly shifted by the 
addition of the short-range attraction. 

A second interesting feature concerns the t dependence of go(r) in the first shell. As the adhesion strength increases 
from t = (HS) to t = 0.8, the contact value monotonically decreases, whereas a discontinuous peak progressively 
builds up at r = (2a)~ . This somewhat counter- intuitive result can be easily understood in terms of the reduction 
of the pressure exerted on particles 1 and 2 by the surrounding ones in the presence of increasing attraction, thus 
providing an average larger separation among 1 and 2. 

Suppose now that we modulate this attraction with the anisotropic dipolar-like dependence described above. When 
ot = 1/2 the effect on g is shown in Fig. [ljb) for three representative values of 6*i 2 : 9yi = (parallel orientation), 
#12 = 7r/2 (orthogonal orientation) and 6*i 2 = 7r (antiparallel orientation). Note that in the orthogonal case the dipolar 
dependence vanishes and one recovers the isotropic behavior. The main differences occur in the first shell, where the 
orthogonal curve (9i 2 = 7r/2 is bracketed between the antiparallel (#12 = 71") and the parallel (6*i 2 = 0) results. 

Similar qualitative results (with different separations among parallel and antiparallel curves) arc found when the 
angle between ?i 2 and Ui is varied. 

From Fig. [TJb) we note that at contact (r = a + ) the antiparallel configuration is more probable that the nose- 
to-tail parallel one; conversely, at separations close to r = 2a~ the parallel alignment is predominant. This can also 
be confirmed by plotting the projections h A (r) and ho{r) of the molecular correlation function /i PY_OL (l,2) on the 
angular basis A(12) and -D(12) respectively. This is depicted in Fig. [2]where the isotropic corresponding contribution 
ho(r) is also reported by contrast. One observes a weak negative correlation for both quantities in the region r w <r + 
and, conversely, a positive correlation close to 2a~ . A crossing occurs approximately around the same value r w 1.7a 
where the parallel component in Fig. [T^b) overtakes the antiparallel one, as expected. 

As we shall see, however, this is a local ordering which docs not affect the condensation process. 

In order to get more insight into such an orientational ordering, we compute the number of particles with orientation 
u 2 that a generic reference particle 1 with orientation Ui 'sees' in an appropriate surrounding volume Vab- Assuming 
Ui as polar axis and taking into account the sphere S with center Ui and radius R, Vab is defined as the portion of 
S corresponding to the solid angle SIab = {(0, f) \0a < < Ob, < ip < 2ir} (see Fig. [3]). Taking for instance 9 a = 
and 9b = 7r/3, we can analyze the 'forward ordering' as seen by the reference particle, while choice 9a = 7r/3 and 
9b = 7r/2 allows to discuss the 'lateral ordering'. 
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The number of particles in an infinitesimal spherical cone of height R = Act and infinitesimal solid angle dfl r in a 
given direction r is dAf (ui, U2, r) = dQ r dr r 2 pg(l, 2), where dfl r = dr. In a finite solid angle Q ab 



7V(ui,u 2 )= / dr f dr r 2 pg(l, 

JQ.AB JO 



2) (13) 



Using the first line of equation ((SJ) and equation (fTTj) we see that, within the PY-OL closure, (g(12)) u Ua = go(r), 
so that the average number is 

JT= (7V(ui,u 2 )) Ui U2 =p 4B / dr r 2 </ (r) - p ^ab <J 3 I (14) 

Jo 

with Qab = Jq 27T Jg^ d9 sin# = 27r(cos# J 4 — cos 6b) and 

7 = (A 3 - 1) /3 + / cfe x 2 /ioacg(x) + A (15) 



Here we have used the results of Rcf. I (see especially Section III D and E), where ho(r) is decomposed into a 
'regular' term ha^ eg (r) and a 'singular' term proportional to the delta function. A similar decomposition is carried 
out (see again in Ref. I) for the fiA(r) and ho(r) parts. Using 

(SIab)- 1 f dr D(l, 2) = M AB (ui • u 2 ) (16) 



Mab = cos 2 9a + cos 9a cos 9b + cos 2 6b — 1 , (17) 

we find that the fraction A" of particles with orientation u 2 in the volume Vab around a reference particle having 
orientation ux, only depends upon the angle 612 = cos _1 (ui • u 2 ) and is given by 

via \ A/"(ui,u 2 ) Ia+MabId , . no . 

X{9 12 ) = = = 1H (ui-u 2 ) (18) 

Al ^0 



/a = J X A ^ x 2 hA, Icg (x) + Aa 
7 D = x 2 h D ^ cg (x) + A D 



(19) 



Fig. 5] (a) depicts X as a function of #i 2 in the case A = 2 (first shell). In the forward region, represented by the 
solid angle &Iab(0, tt/3), we find A"(0) > X(ir) so there are more particles with parallel orientation, with respect to 
particle 1. On the contrary, A"(0) < X(ir) in the surrounding lateral region, characterized by f2 J 4s(7r/3, 7r/2), means 
that here the molecules with antiparallcl orientations prevail. Although these effects are rather small, it is reasonable 
to expect that such differences should grow significantly if the anisotropy parameter at could become much larger than 
the strength t of isotropic adhesion. Note that, while Xforward is larger than ATi a t C rai in the interval < 612 < tt/2, an 
inversion occurs in the region tt/2 < 9 12 < tt in agreement with the results of g(r, #i 2 ) reported above (FigQ] (b)). 

The above results are suggestive of the following physical picture. Because of the limits imposed on the anisotropy 
parameter (0 < a < 1/2) by the choice of the potential, the contribution of the dipolar-like interaction is significantly 
weaker compared to the isotropic part, and does not affect the main condensation process with the formation of 
globule clusters of nearly isotropic shape. This is in sharp contrast with the purely long-range dipolar models which 
are mainly characterized by chain-like aggregatio n 13 i 14 i 15 . However a local ordering occurs within these globular 
agglomerates of condensation, that are mainly formed by short portions of antiparallcl chains running next each other 
and held together essentially by the isotropic attraction. This is schematically depicted in Fig. [4|b). 

We note that only particles belonging to different, antiparallcl, chains have direct contact. Consecutive molecules 
with parallel noisc-to-tail orientation - i.e. belonging to the same chain - are not in contact, but lie with average 
separations slightly smaller than 2ct as suggested by the behavior of g in FigQJb). Thus antiparallel molecules of 
adjacent chains 'mediate' an indirect contact between consecutive particles of a given chain. 

Once again, we stress that this phenomena should be considered a local fluctuation with very short range (of the 
order of one shell, as remarked) and does not extend to the entire fluid. This can be readily checked by considering 
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the limit A — > oo, in which case one finds Ijj = 0, so that the dependence from Ui 112 = cos 6*12 is averaged to zero. As 
we shall see below, a direct consequence of this is that the coexistence line of the isotropic model is not significantly 
affected by the anisotropic part, within the PY-OL approximation. 

In view of the last remark, one might rightfully wonder whether the anisotropic part plays any role in the thermo- 
dynamics of our model. We can convince ourselves that the answer is positive, by considering the exact third virial 
coefficient as defined by 

B3 = ~ j dri dr 2 dr 3 (f (1, 2) / (1, 3) / (2, 3)) Ui , U2 , U3 (20) 

Note that, in view of Eq. ([5|), the exact second virial coefficient B2 = — ^ J dr (/(1,2)) U9 coincides with its 
isotropic counterpart. However, this is not the case for B3, that can be computed following the method outlined in 
Ref. [a for patchy sticky hard spheres, a close relative to the present model. One finds 



where vq = (ir/6) a 3 and 



b 3 = B 3 /v% = 10 - 60t X i + 144t 2 X 2 - 96i 3 X3 , (21) 
Xi = (e(l, 2)} UljU2 X2 = (e(l, 2)e(l, 3)} Ul . U2 . U3 

(22) 

X 3 = (e(l,2)e(l,3)e(2,3)) Ul , U2 , U 3 . 
Again using j5]), we find xi = 1 = X2- The exact value of X3 turns out to be 

11 3 

Xs = 1 - . (23) 

The anisotropic contribution is represented by the term — (ll/72)a 3 s» 0.02, that is very weak with respect to the 
isotropic one. 

Having assessed the limits of the model, we now turn to discuss the limits of the approximation involved in the 
PY-OL closure. A simple and direct way to quantify its deviation from the exact results is to consider the first order 
density expansion of the exact direct corr elation function c(l,2) = /(l, 2) + c^(l, 2)p H .We find 

c«(l,2) = C «_ OL (l,2) + C «(l ) 2) (24) 

with c«_ OL (l,2) = 4 1) (r)+4 1) (r) A(l, 2) + cg } (r) 0(1,2), and 

c«(l,2) = (ai) [ 7 iV) A(l, 2)0(1, 2) +7^ (a) ^(1,2)] aS(r-a). (25) 

Whereas the PY closure includes both Cpy_ OL (l,2) and ci x (1,2), thus reproducing the exact third virial coefficient 
through B 3 = — | J dr {c^ (1, 2)) U1U2 , the PY-OL approximation omits the contributions included in Ccx(l,2). 
Consequently, b^ Y ~ reduces to the purely isotropic contribution: £> 3 so = 10 — 60i + 144t 2 — 96t 3 . The anisotropic 
contribution &| mso ~ ^(at) 3 , stemming from Cox(12), can be easily computed again with the help of Eqs. ([5]), in 
agreement with the exact result Eq. ([23|) . 

Next we consider the thermodynamics. As in all approximate closures, even within the PY one there exist three 
standard routes to the equation of state: compressibility, energy, and virial routes. 

In the first two cases, it is easy to convince oneself that the result is the same as for the isotropic SHS system 
calculated in Ref. [lj|. This is again due to Eq. ([5]) and is a consequence of the linearity of the expansion in the 
angular part involved in the PY-OL approximation, Eq. ([7]), and of the minor role played by the anisotropic part, 
as testified by the weak a— dependence of the third virial coefficient (Eqs. (|2"2"|) and (|2"3"1) ). This is also in agreement 
with the stability analysis of Ref. I, which can also be extended to finite values of the wave vector k. 

As often the case, the virial route is more delicate. Here, standard steps lead to 

y = 1 + 4wb(<0 - 4rft j 22/o (a) + Y '(a) + y[2 to (a) + Y D (a)]} (26) 

where 77 = pvo is the packing fraction, and Y (r) = ry (r), Yjj(r) = ry£>(r), with y PY (l, 2) = 1 + 7(1, 2) being the PY 
cavity function. 
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For given t, r), one can calculate yo(c), ZAd(c), Yq(c) and Y^{a) analytically using expressions from Ref. I. However, 
Yq{o) and Y^(u) require some care, since space derivative and sticky limit do not commute^. So from Eq. (|26p one 
finds the virial pressure. The corresponding results are collected in Fig. [5] for different values of t, at both a = 
(isotropic case) and a = 1/2 (with the anisotropic contribution included). A comparison with the virial expansion up 
to the third virial coefficient is also added in the case t = 0.9 and a = 1/2. 

In agreement with the previous structural findings, we find a dependence on the anisotropy. This is very small for 
t < 0.5 but increasingly appreciable for larger values of the adhesion strength t. In Fig. [S] one clearly sees that the 
pression increases by roughly 10% on going from a = (no anisotropy) to a = 0.5 (maximum anisotropy) for t = 0.9 
and 77 > 0.2. 

Despite the strong differences with the pure dipolar case, it proves instructive to get some insight into the compe- 
tition between the tendency to condensation on the one hand and to chaining on the other hand, by applying to the 
present model the arguments put forward by Tlusty and Safrani^ in the dipolar case. These authors devised a phe- 
nomenological theory, where the two above-mentioned tendencies are represented by the concentrations of 'junctions' 
and 'ends', respectively (see Ref. [13). We have closely followed their arguments to derive the critical parameters r c 
and r\ c of the present model in terms of the energies ex, £3 of ends and junctions respectively. One finds^ 

£i-3e 3 ex (2 In 3 - In 2) -e 3 In 2 

ln?7 c = . (27) 



31n3-21n2 ' ' ex - 3e 3 

which coincides with the results of Ref. Matching this critical values with the one of the isotropic adhesive spheres 
of Miller and Frenkel 2 ^, r c = 0.1133 and r/ c = 0.266, we find ex = 0.186 and e 3 = —0.0102 (the value e 3 < means that 
junctions are enhanced with respect to ends, once more favoring condensation). In our results the number of ends p\ 
and junctions p 3 turn out to be equal (px = P3) at the critical point, which is thus a point of connectivity transition 
in the system. Fig. [6] depicts the coexistence line, which does not display the re-entrance characteristic of the pure 
dipolar case (compare with Fig. 2 of Ref. fl3f ). This is in complete agreement with the remark by Tlusty and Safran 
that the addition of an isotropic short-range attraction - such as the case of the present model - reports the curve 
to its characteristic parabolic shape (see also Fig. 3 of Ref. Il3h . This is also consistent with very recent numerical 
simulation result a 25 ' 26 showing that the addition of a very weak isotropic attraction to the dipolar HS potential makes 
the condensation transition easily observable. 

In summary, we have studied structural and thcrmophysical properties of a particular hard-core fluid where the 
attractive part of the potential includes an anisotropy of dipolar form infinitcsimally short and infinitely strong. 

Any two molecules of the fluid interact only at contact with a potential having, in addition to an adhesive isotropic 
part of the Baxter type, an additional adhesive term, whose intensity depends upon the mutual orientation of the two 
particles in a dipolar fashion. 

Our potential belongs to a class of simple anisotropic models that have recently attracted considerable interest in 
connection with aggregation phenomena in colloidal fluids, polymers and globular proteins, because of their possible 
experimental relevance for self-assembling materials and biological viruses. 

The extremely short-range nature of this peculiar dipolar interaction strongly contrasts with the long-range nature 
of the dipolar hard sphere model. In the latter case, the formation of chains and long anisotropic agglomerates 
significantly affects the possibility of a gas-liquid transition. Using a simplified treatment of the angular part, based 
upon a first-order expansion in angular invariants so to allow an almost fully analytical solution, we have shown 
that only the local (first few) coordination shells are affected by the anisotropy. This is due to the fact that the the 
oricntationally dependent part of the potential has a relatively weak strength with respect to the isotropic attractive 
term, as forced by the particular choice of the potential associated with the < a < 1/2 limits. As a result, all 
structural and thermodynamical properties are only mildly affected by the anisotropic adhesion. 

Nonetheless, the competition of the two adhesive terms (the isotropic and the anisotropic ones) gives rise to an 
anisotropic local ordering within each (almost isotropic) molecular agglomerate consisting of short chains of molecules 
with parallel head-to-tail orientation, 'glued' to similar chains globally oriented in the opposite direction, thus giving 
an antiparallel alignment for particles belonging to two adjacent chains. 

It would be interesting to contrast the present results with more realistic models incorporating a competition between 
an isotropic and anisotropic short range interactions, such as for instance Stockmayer fluids^, dipolar Yukawa HS 
fluids^ or combination of dipolar and square-well potentials^. 

In spite of its simplicity, the results of the present work suggest that, in the presence of dipolar-like anisotropy, 
one can continuously tune from situations only affecting the local ordering (such as in the case presented here) to 
situations where this effect is much more global (such as the real dipolar case), by simply adjusting the range of 
interaction. 
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FIG. 1: (Color online) (a) Isotropic part of the pair correlation function, go{r) = g lsoSHS (r), at rj = 0.4, for t = 0, 0.1, 0.3, 0.5, 0.8 
corresponding to increasing adhesion strength or decreasing temperature, t = yields the HS limit, (b) Behavior of g(r, 612), 
when a = 1/2, at 77 — 0.4 and t — 0.8, for three representative orientations 612 = 0,7r/2, tt (parallel, orthogonal and anti-parallel 
configurations). 




. 2: (Color online) The three components ho(r), /ia(V) and /iu(r) of the molecular total correlation function, for 
0.4, and t = 0.8. At r < a one has ho(r) = —1, /ia(t) = /io(r) = 0. 
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FIG. 3: (Color online) Schematic representation of the volume Vab (included in the shading area) in two different situations: 
(top panel) 6a = giving the contribution from the forward region, and (bottom panel ) 6a 7^ Ob 7^ giving information on 
the lateral adjacent region. 
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FIG. 4: (Color online) (a) Plot of the fraction X of molecules with orientation #12 contained in the volume Vab defined by 
a < r < 2a and by the solid angle Qab, with 8a = 0,8b = 7r/3 for the forward direction and 8a = tt/3,8b = tt/2 for the 
lateral direction. Parameters are: 01 — 1/2, 77 = 0.4 and t = 0.8 in all cases, (b) Schematic representation of a globular cluster, 
with internal chain-like orientational ordering. 
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FIG. 5: (Color online) A plot of fipvo versus the packing fraction r\ for four different values of t = l/(12r), with (a = 1/2) 
and without (a = 0) the anisotropic contribution, using the virial route to the pressure. The 3rd order virial expansion is also 
added (v.e.) in the most relevan case t = 0.9 and a — 0.5 for comparison. The part of the lines which are not shown correspond 
to a loss of solutions. 
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FIG. 6: (Color online) Phase diagram of our dipolar-like sticky hard spheres, calculated for defect energies e\ — 0.186 and 
£3 = —0.0102. At the critical point the binodal curve (continuous line) and the connectivity transition (dashed curve) meet. 
The line denote the end-rich "gas" with the junction-rich "liquid" . 



